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(54) Abstract Htle 

Processing seismic data 

(57) A method of processing seismic data comprising the steps acquiring seismic data in at least two shots 
that have different source-receh/er azimuths from one another, and selecting events arising from interface 
seismic waves in each shot. The interface wave events are processed to provide information about azimuthal 
dependence of the velocity of seismic energy in the near-surface of the earth's interior. 

The method may be used to process events arising from Scholte waves in marine seismic data, or to 
processing events arising from Love waves or Rayleigh waves in land seismic data. 

The method may be applied to seismic data acquired specifically to provide information about the 
near-surface or to seismic data acquired to provide information about the underlying geological structure. 
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FIG 1 





FIG 2 

Circle shoot survey design for azimuthal coverage with the receiver array 
( triangles ) in the centre and the shots ( stars ) located on two circle segments 




FIG 3 

Circle shoot survey design that fills in the gaps in the azimuthal coverage of 
the survey design shown in fig 2. The receiver array is indicated by triangles, 
the shots are indicated by stars and are located on two circle segments 




FIG 4 

Cross-spread survey design for azimuthal dependant phase-velocity 
measurements recorded by a receiver array ( triangles ) . The shots are located 

on lines SI and S3 
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Cross-spread shot pattern for azimathal dependant phase-velocity 
measurements within a source anay ( black stars ).The receivers 
( triangles ) are located on lines Rl and R3 
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Select data in offset range 500m-2000m for receiver 
array, N shots in azimuth range 
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Have all shots been processed ? j 14 
Fit azimuthal terms using equation (4)1- 
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A method of and apparatus for processing seismic data 

The present invention relates to a method of or an apparatus for processing seismic data, 
and in particular relates to a method of or an apparatus for processing seismic data 
acquired at a survey location where the earth's interior at or near the earth's surface has 
different seismic properties to the underlying layer(s) of the earth's interior so as to 
obtain information about the properties of the near-surface region of the earth's interior, 
nie invention also relates to a method of seismic surveying that includes processing 
acquired seismic data in the above vray. The invention may be applied to processing 
seismic data acquired in a marine seismic survey, to processing seismic data acquired in 
a land seismic survey (including data acquired in a permafiost region), or to processing 
seismic data acquired in a bordiole sdsmic survey. 

Figure 1 is a schematic view of one seismic surveying arrangement. This figure shows 
a marine surveying arrangement in which seismic energy is emitted by a seismic source 
1 that is su^ended from a towing vessel 2 beneath the surface of a water column 7, 
such as the sea or a lake. When seismic energy that has been emitted downwards 
reaches the sea-bed. some of the energy will pass into the earth's interior and propagate 
downwards through the earth's interior until it is reflected by a geological feature that 
acts as a reflector 3. The reflected seismic energy passes upwanis through the earth's 
interior, and is detected by an array of seismic receivers 4 disposed on the sea-bed. 
Infomiation about the earth's interior is obtained, for instance, by deteraiining the travel 
time of seismic energy from the seismic source 1 to the receivers 4. Only one reflector 
3 is shown in Figure 1 but. in practice, a number of geological features within the earth's 
interior wiU act as partial reflectors for seismic energy. Moreover, only one source 1 is 
shown in Figure 1 but in practice a marine seismic surveying arrangement may have an 
array of sources. 

The geological structure of the earth is not uniform. One problem encountered in 
processing marine seismic data is that the seismic properties ofthe earth in a region 5 at 
the earth's surface (hereinafter referred to as the "near-surface") are frequently 



or near 
si] 



gnificantiy different from the properties of the underlying geological stracture 6 



(hereinafter referred to as the "basement"). This can occur if, for example, the near- 
surface 5 is less consolidated than the basement. In particular, the yelocity of seismic 
energy may be significantly lower in the near-surface 5 than in the basement 6, and the 
near-surface is therefore sometimes alternatively referred to as a "low velocity layer" 
(LVL). This difference in velocity will produce a shift in the travel time of seismic 
energy compared to the travel time that would be recorded if the near-siuface and the 
basement had identical seismic properties, and these shifts in travel time are generally 
known as "static shifts", or just "statics". 

The near surface 5 is shown at the surface of the earth's intmor in Figure 1(a), but it 
need not extend to the sea-bed and there could be a fiirther layer overlying the near- 
surface. The boundary between the near-surface and the basement inay not be distinct, 
and the earth's properties may change gradually over a finite depthTange. 

The static shift generated by the near-surface 5 depends on the thickness of the near- 
siu-face, and on the velocity of propagation of seismic energy through the near-surface. 
Lateral variations usually occur in both the thickness of the near-surface 5 and the 
propagation velocity through the near-surface, so that the static shift observed at a 
seismic receiver at one location is likely to be different from the static shift observed at 
a receiver at another location. To a first approximation, the entire data set recorded at 
one receiver will be advanced or delayed by a static time shift relative to data recorded 
at another receiver. 

It is highly desirable to take account of the static shift when processing seismic data. 
Unless these static shifts are removed from the seismic data, ambiguity will exist as to 
whether variations m arrival times of seismic events from deeper layers are due to 
variations in the depth or lateral locations of those deeper layers, or simply arise owing 
to propagation effects in the near-surface S. 

The downwardly propagating seismic energy incident on the sea-bed will contain p- 
waves (pressure waves). As the p-waves propagate through the earth's interior it is 
possible for partial mode conversion to an s-wave (shear wave) to occur. In Figure 1, 



partial mode conversion is shown occurring when the seismic energy undergoes 
reflection at a target reflector 3 disposed under the near-surface. (The reflection itself is 
only a partial reflection, since some seismic energy will be transmitted through the 
target reflector without reflection.) As a result the upwardly propagating seismic energy 
produced upon reflection at the target reflector 3 includes both an upwardly prq>agating 
p-wave (shown as a soUd line) and an upwardly propagating s-wave (shown as a broken 
line). 

A wave that undergoes mode conversion ftom a pressure wave to a shear wave at some 
point between the source and receiver is generally known as a '♦ps-wave". The static for 
a ps-wave can be split into a source-side static for the down-going p-wave, and a 
receiver-side static for the upgoing s-wave. 

A further problem in analysing the seismic data is that the seismic properties of the 
near-surface 5 may be azimuthally anisotropic. In this case, an up-going s-wave is split 
as it passes through the near-surface 5 to give a fast s-wave and a slow s-wave which 
have different velocities and particle motions from one another. }b the presence of 
azimuthal anisotropy the seismic data acquired at the receiver 4 mil thus include a pair 
of ps-events. The fast s-wave and the slow s-wave have different polarisation states. 
Furthermore, the anisotropy of the seismic properties of the near-surface 5 means that 
the fast s-wave will have a different static to the slow s-wave. 

R Bale et al have described, in "Analysis and stacking of 3-D converted waves data in 
the presence of azimuthal anisotropy" SEG Expanded Abstracts of 70* Annual Meeting 
(2000), a mefliod of making polarisation measurements on ps-seismic events. After the 
polMisation angle, which is the azimuth of the fast s-wave, has been determined the 
seismic data may be rotated into the fast and slow s-wave directions. 

A Gxst aspect of the present invention provides a method of processing seismic data, the 
method comprising the steps of: selecting events arising from interface seismic waves in 
at least two shots, the two shots having different source-receiver azimuths from one 
another; and processing the selected events thereby to detennine information about 



azimuthal dependence of the velocity of seismic energy in the near-surface of the 
earth's interior. , 

The present invention is thus able to provide information about the near-surface. In 
particular, the invention can provide information about the local phase velocity at a 
point in the near-surface located generally below a receiver. 

The seismic data may be marine seismic data and the selected events may be events 
arising from Scholte waves. 

A second aspect of the invention provides a method of seismic surveying comprising 
the steps of: acquiring first seismic data at a survey location; and processing the first 
seismic data according a method as defined above thereby to obtain information about 
the properties of the near-surface at the survey location. 

The method may further comprise processing the first seismic data to obtain 
information about at least one underlying layer. Thus, it is possible to obtain data 
relating to the basement and data relating to the near-surface iri a single survey* 

Alternatively, the method may fiirther comprise the steps of: acquiring second seismic 
data at the survey location; and processing the second seismic data to obtain information 
about at least one underlying layer at the survey location. This allows the first seismic 
data to be acquired specifically to provide information about the near-surface and allows 
the second seismic data to be acquired specifically to provide information about the 
basement. 

A third aspect of the present invention provides an apparatus for processing seismic 
data, the apparatus comprising: 

means for selecting events arising from interface seismic waves in at least two 
shots, the two shots having different soxirce-receiver azimuths &om one another; and 

means for processing the selected events thereby to determine information about 



azimuthal dependence of the phase velocity of seismic ener^ in the near-surface of the 
earth's interior. 

The jq)paratus may comprise a programmable data processor. 

A fourth aspect of the invention provides a storage medium containing a program for a 
data processor of an ^paratus as defined above. i 

Preferred embodiments of the present invention will now be described by way of 
iUustrative example with reference to the accompanying Figures in which: 

Figure 1 is a schematic iltastration of a marine seismic survey; 

Figure 2 is a plan view of a seismic surveying arrangement suitable for use with the 
present invention; 

Figure 3 is a plan view of a further seismic surveying arrangement suitable for use with 
the present invention; 

Figure 4 is a plan view of a further seismic surveying arrangement suitable for use with 
the present invention; 

Figure 5 is a plan view of a further seismic surveying arrangement suitable for use with 
the present invention; 

Figure 6 is a block flow diagram illustrating a method according^ to one embodiment of 
the present invmtion; 

Figures 7(a) and 7(b) are schematic iUustrations of processmg seismic data using 
information about seismic properties obtained by the present invention; and 

Figure 8 is a block schematic diagram of an apparatus according to the invention. 



The present invention makes use of surface waves or interface waves to obtain 
information about the seismic properties of the near surface of the earth's interior. 



Interface waves are waves that travel in a defined mode along an interface at which a 
velocity contrast exists. One example of an interface wave i^ a Scholte wave, which 
travels along a fluid-solid interface, such as the interface between the sea and the sea- 
bed. A Scholte wave is dispersive, so that its phase- velocity and its group- velocity each 
vary with frequency and with the mode of the wave. As described by R. D. StoU et al, 
in "Measuring Lateral Variability of Sediment Geoaccoustic Properties", J. Acoust 
Soc. Am. Vol. 96, pp 427-438 (1994), a Scholte wave can be excited by a source such 
as an explosive or a shear-wave source disposed at or in the sea-bed. A seismic source, 
such as an air gun, disposed within the water layer can also excite Scholte waves as a 
result of the conversion of seismic energy at the fluid-solid interface. Thus, in the 
seismic surveying arrangement shown in Figure 1 a Scholte wave may be generated in 
the region 8 where seismic energy propagating downwards from the source is incident 
on the sea-bed. A Scholte wave generated by conversion of seismic energy at the region 
8 will then propagate through the sea-bed, close to the interface between the sea-bed 
and the water layer. That is, the Scholte wave propagates through the near surface of 
the earth's interior. ; 

The velocity of a Scholte wave depends mainly on the s-wave velocity of the near 
surface 5, and also on other properties of the near surface such as its p-wave velocity 
and its density. Variations in the s-wave velocity of the near surface 5 thus lead to 
variations in the phase-velocity of a Scholte wave propagating through the near surface. 
Thus, if the Scholte wave phase-velocity can be determined, it can be inverted to obtain 
information about natural variations of the s-wave velocity of the near surface. 

As described by M. L. Smith et al in 'The Azimuthal Dependency of Love and 
Rayleigh Wave Propagation In Slightly Anisotropic Medium", J. Geophys. Res. Vol. 
78, pp332 1-3333 (1973), the azimuthal variation of the velocity of an interface wave in 
the presence of slight azimuthal anisotropy can be described by the following equation: 



+ A,{^)cosM+ A,{a?)smi2f>)+ Aj{a))cos{4^)+ A,{a))sm{4^) (1) 

i 

In equation (1), r(<»,^)is the phase-velocity of the interface waVe at a fiequency © and 
an azimuth Vo(g),(Z') is the reference phase-velocity of the interfece wave as 
calcukited for a standard model. A© (<») is the isotropic component of the difference 
between the actual phase velocity and the value of the phase velocity given by the 
standard model, and Ai to A4 are frequency-dependent quantities that describe the 
azimuthal dependence of the difference between the actual phase velocity and the value 
ofthe phase velocity given by the standard model. • 

The azimuthal direction ofthe surface wave is given by: • 



Pja,, =arctanl — 



(2) 



Figure 1 illustrates the invention as applied to a typical seabed seismic survey. A 
seismic source 1 is towed through the water layer by a towing vessel 2. Although only 
one source is indicated in Figure 1, in practice an anay of separate sdsmic sources, such 
as an array of airguns is generally used. Alternatively, seismic sources such as 
explosives or shear-wave sources, disposed on the sea-bed may l?e used. 

An array of seismic receivers 4 is disposed on the sea-bed. The feceiver array may 
conveniently be embodied as an ocean bottom cable deployed on the sea-bed by another 
survey vessel 2' but, in principle, the seismic receivers 4 could be deployed individually 
on the sea-bed. Each seismic receiver may be a single-component seismic receiver, that 
records only a single component ofthe seismic wave-field at the receiver. Each seismic 
receiver 4 may alternatively be a multi-component seismic receiver ttiat can record more 
than one component of the received wave-field. In particular, each seismic receiver 4 
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may include a three-component geophone for recording three orthogonal components of 
the particle motion (the term "particle motion" includes particle displacement, particle 
velocity, particle acceleration, and higher derivatives of the particle velocity), and may 
possibly also include a hydrophone for recording the water pressure at the sea-bed 
(which is a scalar quantity). 

I 

In principle the invention may be applied to determining the static shift of the near- 
surface from a full seismic survey - that is, from a seismic survey that is intended to 
provide information about the basement 6. However, a seismic survey intended to 
provide good quality data relating to the basement 6 may not always provide good 
quality data relating to the near-surface 5. Therefore, a separate survey to determine the 
static shift of the near surface can be carried out before or after a fiiU seismic survey is 
carried, since this might allow the static shift to be determined more accurately than if 
the static shift is determined from the full survey. The present invention can be applied 
to seismic data acquired in a full seismic survey, or to seismic data acquired in a seismic 
survey cairied out specifically to provide information about the static shift of the near 

surface. i 

I 

Scholte waves may be generated in the seismic survey of Figure 1 , 'as a result of 
conversion of oieigy at the fluid-solid intrafece (i.e. at the interface between the water 
layer and the near sur£u;e 5). It is assumed that the Scholte waves are generated in the 
region 8 generally below the source array 1. 

When the seisanic source array 1 is actuated, the seismic data acquired at each of the 
receivers 4 will contain events relating to at least the following: 

(a) seismic energy that has travelled direct fiom the source array 1 to the receiver 4 
through the water column (the "direct wave", not shown in Figure 1); 

/ 

(b) p-wave arrivals that have undergone reflection at the reflector 3 below the near 
siu-face 5; 
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(c) s-wave arrivals arising fiom reflection and partial mode conversion at the 
reflector 3 below the near surface 5; 



(d) intCTface waves such as Scholte waves; and i 

(e) events arising from reflection, and possible mode conversion, at other reflectors 
located within the basement 6 (the method of the invention does not make use of these 
events, but they will generally be present in the seismic data). 

It should be noted that the amplitude of Scholte waves depends primarily on the 
distance between the seismic source and the sea-bed. In shallow Water, seismic data 
acquired at the receivers is likely to contain adequate Scholte wave data, but in deep 
water the Scholte wave events may have a low signal to noise ratio. If the seismic data 
acquired at the receivers 4 does not contain Scholte wave data having a good signal to 
noise ratio, it may be necessary to carry out the additional survey described above. The 
additional survey can be specifically focused on the excitation and recording of Scholte 
waves in order to ensure reliable results. 

Once the seismic data has been acquired, it is initially necessary to select one 
component of the seismic data for the purpose of determining the phase-velocity of 
seismic energy in the near-surfece 5. In a fuU 4-C survey, data relating to the vertical 
component and to orthogonal horizontal components of the parti4le motion at the 
receiver, and data relatmg to the pressure in the water colunan wi^l be available. It is 
generally preferable to select the vertical component of the particle motion for further 
analysis, although data relating to the pressure in the water colunin may alternatively be 
selected. 

The horizontal components of the particle motion can suffer from geophone coupling 
problems, which give rise to possible phase and ampUtude variations within the receiver 
array. If the horizontal components of the particle motion are selected it is preferred, 
although not essential, to rotate the co-ordinates towards the azimuth ^of the source. 
This may be done using: • 
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X" 






-sin^ cos^ 





(6) 



where Vx and ^ are the measiired x- and y-components of the particle motion, and Vxg^ 
and Vyp are the components after rotation by q>, \ 

K 

\ 

The number of receivers in the receiver array should be selected depending on the 
required data quality and resolution. Typically, an array of from 4!to 20 receivers has 
been found satisfactory. I 

The shots for which the data is to be analysed are selected to provide as wide a possible 
range of source-receiver azimuth. The shots are selected so that the minimum offset 
from the source array to the centre of the receiver array is several times greater than the 
wavelength of the Scholte waves. The maximum offset used is determined by the 
amplitude decay of the Scholte wave, and by the recording time. For example, a 
Scholte wave propagating at 2Hz typically has a phase- velocity of around 300m/s, 
resulting in a wavelength of around 1 50m. A suitable minimum offset would be 
around 500m, which is slightly over three times greater than the wavelength. In the 
case of a IDs recording length, the maximum offset should be 300Qm. In order to avoid 
spatial aliasing, a maximum geophone separation of around 25m is suitable, although 
this will depend on the minimum expected wavelength of the Scholte wave. 

Once suitable shots have been selected, the phase-velocity of the Scholte wave is 
determined for each shot. In a preferred embodiment, this is done by calculating the 
local frequency-wavenumber spectrum ("fk-spectrum**) for each shot. The peaks in the 
fk-spectrum correspond to the phase velocity of the Scholte wave. ; The local fk- 
spectrum will typically contain more than one peak, and different peaks can represent 
different Scholte wave modes or other signals. The velocity that is detemiined is a local 
measiirement of the velocity in a region below the receiver used to acquire the shot. 



The fk-spectnim may be determined by using, for example, the "MUSIC" (multiple 
signal classification) algorithm developed by H. Krim and M. Viberg in 'Two Decades 
of Array Signalling Research", IEEE Signal Processing Magazine, July 1996 and by 
K. Iranpour and E. Muyzert in "Application of Parametric Methods to a Number 
Estimation in Seismic Data", submitted to Geophysical Prospecting. The "MUSIC" 
algorithm is a parametric method that makes it possible to estimate an fk-spectrum from 
a small number, typically 4 to 20, of equi-distance spaced receivers. The fk-spectrum is 
a local measure, since the dispo^ion is measured within the receiver array and is 
indqjendent of the distance of the source. The phase velocity of jthe Scholte wave is 
givai by the maxima in the fk-spectrum. At each frequency (o the maximum energy in 
the local fk-spectrum is picked manually or automatically. The i>hase-velocity c is 
calculated from the wavenumber corresponding to maximum energy using c = a>lk. 
When more than one distinct peak is found in the fk-spectrum it is necessary to identify 
which peak is which mode. This can be done by comparing the phase-velocity 
estimated from the fk-spectrum with the data in a Jime domain representation. 

In the case of a receiver array in which the receivers are disposed substantially along a 
straight line, and for a shot in which the source-receiver azimuth'lies along the axis of 
the receiver array, the local phase-velocity is equal to the apparent phase-velocity 
obtained from the fk-spectrum. However, if the source-receiver 'azimuth of the shot 
makes an angle a with the axis of the receiver array, it is necessary to correct the 
apparent phase-velocity determined from the fk-spectrum to account for this. In this 
case, the tme phase velocity of the Scholte wave is given by: j 

V^{a^,a)^ Vjj, {a>.a)/cos a (4) 

where Vtmc {<^y ot) is the true phase-velocity of the Scholte wave, and Vfk(co, a) is the 
apparent phase velocity determined from the fk-spectrum. ' 

Equation (1) above may be solved for the coefficients Ao, Ai, . . ;..A4 if the phase- 
velocity of the Scholte wave is known for at least five different azimuths at a particular 
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frequency. In this case. Equation (1) may be solved to give the value of the coefficients 
Ao. . • A4 at the particular frequency. 

\ 

To solve equation (1) it may be re-written as the following matrix equation that is to be 
solved for the imknown coefficients Ao. . . . A4: 



1 cos 2^ I sin 2^ J cos4g>^ sin4f?j 
J cos2f>j^ sm2^j^ cos4<Pfj cos4^j^ 



aXo>) 



(5) 



In equation (5), <p\ . . ..(z;sf represent the different azimuth angles, N rnust be 5 or greater, 
©j represents one particular frequency, and Vi (©j, ^) represents tlje phase-velocity 
determined for the Scholte wave at frequency coj for a shot having a source-receiver 
azimuth . As noted above with respect to equation (1), Vo(a)j, ^) is a reference 
phase-velocity calculated for frequency coj and source-receiver azimuth ^ using a 
standard model. Details on how to derive a suitable reference model are given by E. 
Muyzert, in "Schohe wave inversion for a near surface S-velofcity model and PS-statics" 
SEG Expanded Abstracts ppl 197-1200 (2000). 

Equation (5) may be solved to give the values of the coefficients Ap to A4 at frequency 
coj, provided that the phase-velocity at frequency ©j is known for at^ least five different 
source-receiver azimuths. Equation (5) may be solved by any suitable method, for 
example a least-squares method. 



Once the coefficients Ai (coj) are known, the polarisation angle at frequency coj may be 
determined using equation (2). 



" i 

Once the coefficients Aj have been detennined for a frequency ©i, the process of solving 
equation (5) to detennine the coefficients Ai may be repeated at one or more different 
frequencies. This provides the values of each coefficient Aj at two or more different 
frequencies. Furthermore, the polarisation angle may be determined from equation (2) 
at each frequency for which the value of the coefficients Aj ((d) is known. 

Alternatively, Equation (5) may be solved just for the coefficient^ Ao, A, and A2 since 
the other coefficients are not requued for the calculation of the polarisation angle. This 
is an approximation, but requires measurements at only three or '^ore different azimuths 
(that is, N ^ 3 instead of N ^ 5). Equation (5) may be solved just for the terms Ao, A, 
and A2 using: ■ 



1 cos 2^, sin 2^, 
1 cos 2^^ sin 2^^ 



AM) 



(6) 



nie static shift of the near surface 5 may be estimated by deteraiining the frequency- 
average value of the polarisation angle, using: I 



9>ftai iaverage)= J^ll^M K ) 



(7) 



In equation (7), the average polarisation angle is determined from the values of the 
polarisation angle forM different values of the frequency «>j. Alternatively rather than 
averaging the polarisation angles obtained at different frequencies, tiie polarisation 
angle ^ corresponding to a frequency that samples the near surface well, such as the 
lowest available frequency, could be used. 

J 

The average polarisation angle obtained using equation (7), or tj»e polarisation angle 
corresponding to a frequency that samples the near surface |well. may be used to 



remove the effects of the anisotropy of the near surface 5 from ps-converted wave data. 
This may be done by ""rotating" the acquired data such that one horizontal component is 
coincident with the polarisation angle. The "rotation" is done using an equation similar 
to equation (5) above, but with ^in equation (5) replaced by ^ast or ^^st (average). 

This process is illustrated schematically in Figures 7(a) and 7(b). Figure 7(a) shows the 
orientation of the two horizontal internal axes of the receiver, which are assumed to be 
the receiver's x- and y-axis, and their relation to the frequency-averaged polarisation 
angle. It is assumed that there is an angle p between the receiver's x-axis and the 
average polarisation angle. 

In Figure 7(b) new axes have been defined for the receiver. These 'are shown as the x'- 
and y'-axes. The x'-axis is at an angle (3 to the original x-axis of Figm-e 7(a), and so is 
coincident with the polarisation angle. As a resuh, the x '-component of the receiver 
data now contains the fast ps-wave, and the y' component of the seismic data acquired 
at the receiver contains the slow ps-wave. 

The process of rotating the receiver axes is a virtual rotation process. The x- and y- 
components of the seismic data acquired at the receiver are processed to produce the 
data that would have been acquired if the receiver axes had been aligned along the x'- 
and y '-directions of Figure 7(b). i 

I 

Once the horizontal components of the seismic data have been rotated so that one 
component is coincident with the polarisation angle, the static shift may then be 
detemiined for both the x'-component of the seismic data (that is, for the fast ps*wave 
data) and for the y' component of the data (that is for the slow ps-wave data). This may 
be done using any conventional technique for determining the static shift. 

Once the static shift has been determined for both the fast and slow ps-waves, the arrival 
time of the slow s-wave may be shifted by the difference between the static shift 
determined for the fast ps-wave and the static shift determined for the slow ps-wave. 
The x'- and y' -components of the seismic data may now be rotated back by an angle --IJ, 
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to return to the original x- and y- axes. The effect of the azimttthal anisotropy of the 
low velocity layer has now been removed from the ps-wave data,;since the difference in 
static shift between the fast ps-wave data and the slow ps-wave data has been 

eliminated. ; 

i 

i 

Figure 6 is a block flow diagram showing the principle steps of a method of deteraiining 
the polarisation angle using the present invention. 

At step 10, seismic data having a suitable ofifeet range are selected. The seismic data 
relate to N shots, each shot having a source-receiver azimuth within a desired nmge. 

The present invention may be applied to pre-existing seismic dat^, in which case step 10 
may optionally be preceded by the step (not shown) of retrieving seismic data fiom 
storage. Alternatively, the present invention may be performed on seismic data that 
have been specifically obtained for this purpose, in which case st^) 10 may be preceded 
by the step (not shown) of acquiring seismic data. j 

At step 1 1 the local fiequency-wavenumber spectrum (fk spectrum) for the n* shot is 
determined for the receiver array. This may be done, for example, using the "MUSIC" 
algorithm. 

At step 12 the apparent phase-velocity of the interface wave, in this embodiment a 
Scholte wave, is determined fiom the fk spectrum, as corresponding to the peak value in 
the fk spectrum. It might be possible to determine the apparent i?hase velocity for two 

or mode modes of the interfece wave in this step, if the fk-spectrum has two or more 

I 

maxima.. 

If the n* shot has a source-receiver azimuth that is not aUgned with the axis of the 
receiver array, the apparent phase-velocity is corrected for this at step 13, using 
equation (2). If the source-receiver azimuth of the n* shot is aligned vrith the axis of the 
receiver airay, the apparent phase-velocity obtained in step.l2 will not change as a 
result of the correction in step 13. 
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At step 14 a deteraiination is made as to whether all the N shots selected at step 10 have 
been processed according to steps 1 1, 12 and 13. If step 14 yields a "no" 
determination, another shot is selected and steps 1 1 to 14 are repeated for the new shot. 
Steps 1 1 to 14 are repeated until a "yes" determination is obtained at step 14. 

i 

At step 15 the coefficients A,(a)) in equation (1), or at least the coejfficients Ao, Ai and 
A2, are determined. This may be done by solving equation (5), or by solving the 
reduced equation (6). The result of step IS is, if equation (5) is used, the value of the 
coefficients Aj at a particular frequency ooj or, if equation (6) is usdd, the values of Ao, 
Ai and A2 at a particular firequency ooj. 

At step 16 the polarisation angle at frequency coj is deteraiined using equation (2) and 
the values of Ai((ji)j) and AaCcoj) determined at step 15. 

At step 17 a different frequency is selected, and st^s IS and 16 are then repeated for 

li 

the new frequency. Steps 17, IS and 16 may be repeated for further frequencies. 

*■ i 

At step 18 the mean polarisation angle is determined using equatioh (5). 

Once step 1 8 has been carried out the average polarisation angle is known for the 
receiver array and it is possible to detemfiine the static shift and remove anisotropy from 
the data acquired by that receiver array. In practice, howevCT, it is more efficient to 
determine the average polarisation angle for all receiver arrays before determining the 
static shifts and removing the anisotropy from the data. In the' embodiment of Figure 6, 
therefore, at step 19 a determination is made as to whether all receiver arrays have been 
processed. If step 19 yields a "no*' determination, steps 10 to 18 are repeated for 
another receiver array. Steps 10 to 18 are repeated until a "y^** detjermination is 
obtained at step 19. 



Once a •'y^s" determination is achieved at step 19, it is now possible to use the average 
polarisation angle to remove the effect of anisotropy of the near Surface from the 
seismic data, as described above. 

Figure 2 illustrates one possible survey geometry suitable for acquiring seismic data that 
can be processed according to a method of the present invention. In the survey 
geometry of Figure 2 the receivers 4 are represented by solid triangles, and the sources 1 
are represented by stars. 

In the survey geometry of Figure 2 the receivers 4 are arranged i^ a linear array. Six 
receives are shown in Figure 2, but the survey geometry is not limited to a receiva: 
array having exactly six receivers. The numbor of receivas may be selected to ensure 
that the quality and resolution of the data are satisfactory. In practice, between 4 and 20 
receivers are generally found to give satisfactory results. ; 

In tfie survey geometry of Figure 2 the sources are located generally on the diameter of 
a circle that is centred on the centre of the receiver array. 34 sources 1 are diown in 
Figure 2, but the survey geometry is not limited to Ais exact number of sources. All the 
source-receiver pairs have ^proximately the same offset, ahd this is approximately 
equal to the radius of the circle on which the sources one are arranged. As explained 

above, the oflfeet should be at least a few times the expected wavelength of the Scholte 

1 

wave, and the maximum offeet is determined by the amplitude decay of the Scholte 
wave and by the time period in which data is recorded after a source has been actuated. 

! 

In operation, each source is actuated in turn. When a source is actuated data will be 
acquired at each receiver 4 in the receiver anray, and the apparent phase-velocity of the 
Scholte wave is determined frcrai flie seianic data as described above. The apparent 
velocity is then corrected to take account of the angle between the source-receiver 
azimuth and the axis of the receiver array. 

i 

In principle, the sources can be disposed around the circumference of the circular shot 
line 9. In practice, however, if the angle a between the source-receiver azimuth and the 
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shot line has a magnitude that is greater than around 60^, the apparent phase- velocity is 
so high that the tme phase velocity obtained using equation (2) becomes less reliable. It 
is therefore preferable if the sources 1 are disposed such that the angle between the 
source-receiver azimuth and the axis of the receiver array does not ^significantly exceed 
60^. This has been achieved in Figure 2 by arranging the sources ij so that they have an 
azimuth in the range -60"^ < q>< 60^ or 120"^ ^< 240^. (In Figur6 2 the receiver axis 
is coincident with the 0 axis, so that the source-receiver azimuth of a shot is, 
neglecting the finite extent of the receiver array, equal to the azimuth of the receiver.) 

The seismic surveying arrangement shown in Figure 2 will produce seismic data having 
two gaps in its azimuthal coverage, since there will be no data for azimuthal angles 
between 60'' and 120° or between 240"^ and 300**. To overcome this, it is possible to 
combine the seismic data acquired in the seismic surveying arrangement of Figure 2 
with seismic data acquired in a complimentary seismic surveying arrangement shown in 
Figure 3. This again contains a linear array of seismic receivers 4, but the axis of the 
receiver an^y is aligned along the azimuth ^= 90*". The seismic sdurces are again 
arranged on the circumference of a circle that is centred on the centre of the receiver 

array. The sources are again disposed such that the angle a between the source-receiver 

i 

azimuth and the axis of the receiver array does not exceed 60**, This is achieved by 
disposing the sources 1 in the azimuth range 30° < ^< 150° and 210° < ^< -30°. As in 
the seismic surveying arrangement of Figure 2, the seismic surveying arrangement of 
Figure 3 is used to acquire data by actuating each source 1 in turn. When a source 1 is 
actuated it emits seismic energy and this is recorded by the recbivers 4 of the receiver 
array, and the phase-velocity of the Scholte wave is determined as described above. 

![ 
I 

It will be noted that the seismic surveying arrangement of Figure 2 provides information 
on the phase velocity of Scholte waves for azimuths in the range -60° < ^ < 60° and 
120° < ^< 240°, whereas the seismic surveying arrangement of Fi^e 3 provides 
information about the phase-velocity of Scholte waves in the azimuth ranges 30° £ q>< 
150° and 210° < -30°. Thus, the two seismic surveying arrangements together 
provide coverage of the entire azimuthal range. Furthermore, the azimuthal coverage of 
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the seismic surveying arrangement of Figure 2 overlaps with the pmuthal coverage of 
the seismic surveying arrangement of Figure 3, for example in th? range 30° < ^< 60°. 
These regions of overlapping azimuthal coverage may be used to determine that the 
results of the two seismic surveying arrangements are consistent with one another, and 
also to provide an estimate of the errors in the measurement of the phase velocity. 

It would in practice be possible to embody the seismic surveying arrangements of 
Figures 2 and 3 in a single seismic survey in which receivers were disposed around a 
circular shot line 9 so as give complete azimuthal coverage. For example, a source one 
could be disposed every 10° around a circular shot line 9. the seismic surveying 
arrangement would have receivers that provided both a linear receiver array along the 
<p = 0" line and could also provide a linear receiver array along the fine ^= 90°. This 
could be achieved by having receivers arranged in a cross-shaped receiver array, with 
the centre of the cross being substantially coincident with the centre of the circular shot 
line 9. Alternatively, the receivers could be arranged ina grid - a6x6 grid of receivers 
would provide the receiver array of Figure 2 and the receiver array of Figure 3. Such a 
grid may conveniently be embodied by neighbouring OBC sensor arrays that are 
disposed parallel to one another and close to one another (for example with a 25m 
separation between adjacent OBC arrays). When a particular sei^ic source is actuated, 
seismic data acquired by the appropriate receivers is selected forjanalysis, so that the 
seismic surveying arrangement of Figure 2 or Figure 3 may be r^roduced. 

Figure 4 shows a further seismic surveying anrangement that may be used to provide 
seismic data that can be processed by the method of the present invention. This is the 
well-known "cross-spread" geometry, which comprises two shot lines 20. 21 which are 
parallel to one another, and which are peipendicular to a receiver line 22. The spacing 
between adjacent sources 1 on a shot line 20, 21 is smaU. for example around 25m. The 
spacing between a shot line 20, 21 and the receiver array is large, and may typically be 
several hundreds of metres, for example 600m. As in the seismic surveying 
arrangements of Figures 2 and 3, it is preferable that the angle bptween the source- 
receiver azimuth and the axis of the receiver array does not sigiuficantly exceed 60°, to 
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ensure reliable results are obtained. This may be achieved by choosing the length of the 

shot lines and their separation from the receiver array appropriately. 

i 

1 
i 

Unlike the seismic surveying arrangement of Figures 2 and 3, the pfTset of the source- 
receiver pairs is not constant in the seismic surveying arrangement' of Figure 4. The 
length of the shot lines 20, 21 and their separation from the receiver array is preferably 
selected to keep the offset between the minimum offset and maximum offset described 
above. 

The seismic surveying arrangement of Figure 4 again has gaps in azimuthal coverage, 
since there is no coverage for azimudis between 60^ and 120*^; or for azimuths between 
240** and -60^. In principle, these gaps in the azimuthal coverage could be overcome by 
using a complementary cross-spread surveying arrangement in whiich the receiver array 
was aligned vertically and the source arrays were aligned horizontally. 

A 

i 

Figure 5 illustrates an alternative seismic surveying arrangement that can be used to fill 
in the gaps in the azimuthal coverage of the seismic surveying arrangement of Figure 4. 
This makes use of the principle of reciprocity, which states that the results of a seismic 
survey are unaffected if the positions of the sources and the positions of the receivers 
are interchanged. 

Thus, the seismic surveying arrangement shown in Figure 5 contains a linear array of 
sources one, and two parallel linear arrays of receivers 4. The receiver arrays 23, 24 are 
parallel to one another and are perpendicular to the axis of the source array 25. The 
centre of the source array 25 is approximately midway between one receiver array 23 
and the other receiver array. The phase velocity of the Scholte wave is measured over 
the source array, so that the local phase-velocity of the Scholte wav'p at a point in the 
near-surface 5 underneath the sources 1 is obtained. This £q)proach: assumes that the 
Scholte wave conversion point is close to the position of the source array, and that there 
are no significant variations in the conversion of p-waves to Scholte waves within the 
area underneath the source array. 
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As is the case for other seismic surveying arrangements, the seismic data produced by 
the seismic surveying arrangement of Figure 5 will enable the apparent phase-velocity 
of the Scholte waves to be determined. In order to obtain the true phase-velocity of the 
Scholte wave it is necessary to correct for the angle between the source-receiver 
azimuth and the axis of the source array, by using equation (2), 

When a seismic surveying array according to Figure 4 and a seismic surveying array 
according to Figure 5 are used, centred on the same point, the local phase-velocity of 
Scholte waves is sampled at the same location but at different azimuths. 



The method of the invention has been described above witii respect to marine seismic 
data, and with respect to using Scholte waves to determine the statics of the near surface 
5. The invention is not, however, limited to use with marine seismic data. It may be 
applied to seismic data acquired in a conventional land-based seismic survey, or to 
borehole seismic data acquired in a seismic survey in which the sources or receivers are 
located in a borehole. In the case of a land-based seismic survey interface waves, 
namely Love waves or Rayleigh waves, are generated in and propagate through the near 
surface, and these interface waves may be used to determine information about the 
properties of the near surface. In the case of a borehole-ba^ed seismic survey interface 
waves, namely Stonely waves, are again generated in and propagate through the near 
surface, and the method of the invention may be applied to these interface waves. 

I 

Figure 8 is a schematic block diagram of an apparatus 26 according to the present 
invention. The apparatus is able to carry out a method according to the present 
invention. 

The apparatus 26 comprises a prograitunable data processor 27 with a program mraiory 
28, for instance in the form of a read only memory ROM, storing a program for 
controlling the data processor 27 to process seismic data by a method of the invention. 
The system further conq)rises non-volatile read/write memory 29 for storing, for 
example, any data which must be retained m the absence of power supply. A •forking" 
or "scratchpad" memory for the data processor is provided by a random access memory 
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(RAM) 30. An input device 3 1 is provided, for instance for receiving user commands 
and data. An output device 32 is provided, for instance for displa>|ing information 
relating to the progress and result of the method. The output devi<ie may be, for 
example, a printer, a visual display unit or an output memory. 



Seismic data for processing may be supplied via the input device 31 or may optionally 
be provided by a machine-readable store 33. 

The program for operating the system and for performing the method described 
hereinbefore is stored in the program memoiy 28, which may be embodied as a semi- 
conductor memory, for instance of the well-known ROM type. Hcjwever, the program 
may be stored in any other suitable storage medium, such as magnetic data carrier 28a 
(such as a "floppy disc") or CD-ROM 28b. 
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CLAIMS: 

1 . A method of processing seismic data, the method comprising the stqjs of: 
selecting events arising from interface seismic waves in at least two shots, the two shots 
having different source-receiver azimuths from one another; and processing the selected 
events thereby to determine information about azimuthal dependence of the velocity of 
seismic energy in the near-surface of the earth's mtenor. j 

2. A method as claimed in claim 1 whoein the seismic data ^are marine seismic 
data and the selected events are events arising from Scholte waves. 

3. A method as claimed in claim 1 whoein the sdsmic data are land seismic data 
and the selected events are events arising from Love waves or Rayleigh waves. 

4. A method as claimed in any preceding claim further comprising the step of 

determining the local phase velocity of the interface wave for each shot. 

> 1 
1 

5. A method as claimed in any preceding claim and comprispg the step of 
determining the local phase velocity of the interfece wave at a point in the near-surface 
located generally below a receiver. 

6. A method as claimed in claim 4 or 5 wherdn the step of determining the phase 
velocity of the interface wave in the surface layer for each shot comprises determining 
the local frequency-wavenumber spectrum (fk spectrum) for each shot. 

7. A method as claimed in claim 5 or 6 and further comprising the step of 
correcting the phase velocity in the near-surface of the interface jvave determined for a 

shot for the source-receiver azimuth of the shot ^ 

i 

8. A method as claimed in claim 7 wherein the seismic data] are seismic data 
acquired using a substantially linear receiver anray, and the step of correcting for the 
source-receiver azimuth of a shot comprises determining: 
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Vcorr i^^d) ^V{a^a)/ cos a 



where a is the angle between the source-receiver azimuth of the shot and the line of the 



receiver array. 

9. A method as claimed in any preceding claim and comprising determining the 
phase velocity of the interface wave at a first frequency in the near-surface for at least 
five different source-receiver azimuths; and determining the value of at least the 
coefficients A© to A2 at the first frequency using j 



where V{ii}, ^)is the phase-velocity of the interface wave at a frequency <o and an 
azimuth ^, Vo(o>,^) is a reference phase-velocity of the interface wave, and Ao (co) is the 
variation in the isotropic phase-velocity. 

10. A method as claimed in claim 9 and further comprising the step of determining 
the azimuth ^ast of a fast interface wave at the first frequency using 



11. A method as claimed in claim 8 and fiirtfaer comprising determining the azunuth 
^ of the fast interface wave at at least a second frequency different from the first 
firequency. 

12. A method as claimed in claim 1 1 and comprising the step of determining the 



+ A,{^)cos{2^)-^ .4^(fi>)5in(2^)+ As{a?)cos{4^)'\' A^iio). 





(2) 



average value of the azimuth ^ast of the fast interface wave at at least the first and 
second frequencies. * j 
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13. A method of seismic surveying comprising the steps of: acquiring first seismic 
data at a survey location; and processing the first seismic data according a method 
defuied in any of claims 1 to 12 thereby to obtain infonnation abcjut the properties of 
the near-surface at the survey location. 'j 

i 

1 

14. A method of seismic surveying as claimed in claim 1 3 and further compnsmg 
processing the first seismic data to obtain information about at least one underlying 
layer. 

15. A method of seismic surveying as claimed in claim 13 and further comprising 
the steps of: acquiring second seismic data at the survey location; and processing the 
second seismic data to obtain infonnation about at least one underlying layer at the 
survey location. j 

16. An apparatus for processing seismic data, the ^paratus cpmpnsmg: 

means for selecting events arising from interface seismic jvaves in at least two 
shots, the two shots having different source-receiver azimuths frcjm one another; and 

means for processing the selected events thereby to determine information about 
azimuthal dependence of the phase velocity of seismic energy in the near-surfece of the 
earft's interior. 

17. An apparatus as claimed in claim 16 and fiuAer comprising means for 
determining the local phase velocity of the interface wave in the near-surface for each 
shot. i 

18. An apparatus as claimed in claim 17 wherein the means for determining the 
local phase velocity of the interfece wave in the near-surface for jcach shot comprises 
means for determining the firequency-wavenumber spectrum (fk spectrum) for each 
shot. 
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19. An apparatus as claimed in claim 17 or 18 and further comprising means for 
correcting the phase velocity of the interface wave in the near-surface detennined for a 
shot for the source-receiver azimuth of the shot. 



20. An apparatus as claimed in claim 19 wherein the means for correcting the phase 
velocity of the interface wave are adapted to determine 



Vcorr i^yCt) = V(a),a) J cos a 



where a is the angle between the source-receiver azimuth of the shot and the line of a 
substantially linear receiver array. 

21. An apparatus as claimed in any of claims 16 to 20 and comprising means for 
determining the phase velocity of the interface wave at a first frequency in the near- 
surface layer for at least five different source-receiver azimuths; and determining the 
values of at least the coefficients Ao to A2 at the first fi^equency using 



i 

+ A, {a?)cos{2f>) + A^ {a})sin{2^) + A^ {{i))cos{4^)'h A^ {a})sm{4^) 

i 
j 

where V{a?, ^)is the phase-velocity of the interface wave at a frequjency © and an 

azimuth ^, Vo(a),^) is a reference phase-velocity of the interface wave, and Ao (©) is the 

I 

variation in the isotropic phase-velocity. j 



22. An apparatus as claimed in claim 21 and further comprising means for 
determining the azimuth of a fast interface wave at the first frequency using 



(P/axi=arctan|^-^j 



(2) 
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23. An ^paratus as claimed in claim 22 and further comprising means determining 

the azimuth ffest of the fast interface wave at at least a second frequency different from 

i 

the first frequency. ^ 

24. An apparatus as claimed in claim 23 and comprising means for determining die 
average value of the azimuth ^ of the fast interface wave at at least the first and 
second fi:«quencies. 

25. An {q>paratus as claimed in any of claims 16 to 24 and comprising a 
prograrrunable data processor. 

26. A storage medium containing a program for a data processor of an apparatus as 
defined in claim 25. | 
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